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A GENERALIZATION OF CUNNINGHAM’S EXTENSION OF 
STOKE’S LAW FOR THE FORCE ON A SPHERE 


GALE YOUNG 
THE UNIVERSITY OF CHICAGO 


Cunningham’s formula for the force on a sphere moving within 
a larger concentric spherical boundary is extended to cover a gen- 
eral state of motion of the fluid between them. 


Stoke’s formula for the force on a sphere moving slowly through 
a viscous fluid of infinite extent has played an important role in many 
physical investigations, and has been the subject of much experimen- 
tal and theoretical study. Various modifications of it have been 
worked out by different authors; for some account of these see Smo- 
luchowski (1913), Lamb (1924, p. 565-583, 608), and Bulletin of the 
National Research Council (1931, Chapter 7). 

Cunnigham (1910) worked out the force on a sphere moving 
through a fluid enclosed within a larger concentric sphere, and from 
this tried to estimate the effect of the presence of many other moving 
spheres. His results have been used (Heilbrunn, 1928, Chapter 5) in 
determining the viscosity of a fluid from observation of the rate of 
motion through it of a large number of particles under the action of 
applied forces or thermal agitation. For a discussion of Cunning- 
ham’s work and the many sphere problem see Smoluchowski (1913). 

Cunningham found the force on the inner sphere under the as- 
sumption that the fluid is at rest* over the surface of the outer 


sphere, and then considered various sizes for the larger sphere to 


somewhat correspond to the presence of other particles. Since there 
is in such a situation no surrounding sphere where the fluid is at rest, 
and since it is quite easy to remove the assumption of an exactly speci- 
fied flow at either sphere, it may be worth while to note the general 
result. This will also be of use in other connections. 


I 
The complete solution of the equations of slow, potential force 
motion of an incompressible fluid is given by Lamb (p. 562-564) ina 
form adapted to a region with concentric spherical boundaries. In 
* Williams (1915) treated the same concentric sphere problem by a different 


method. We shall here more nearly follow Cunningham. Williams gives some 
diagrams showing flow patterns, both theoretical and observed. 
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such a motion the mean pressure p of the fluid is a harmonic function 
and may thus be written as 
P= Pn (1) 
where p, is a solid harmonic of degree n. The summation is over n 
and over the various independent solid harmonics (there are 2n + 1 
of these of positive degree n or of negative degree —n — 1) of each 
degree. The fluid motion cannot give rise to a pressure term of de- 
gree —1, so that p., is always zero. Denoting by u,v, w the compo- 
nent velocities of the fluid along the (x, y, z)-axes respectively we 
have 
os is cP Tete 0 Pn 
no |2(20-+1).0% | (n+1)(2n4+-1)(2n5-8) 0a | 
(2) 
Ogn  O%n 0 Xn 
3 es ihe sateen Lie 
+ Ox “rath oy verse st 
with the corresponding expressions for v and w obtained by cyclical 
interchange of x, y, and z. Here y is the coefficient of viscosity of 
the fiuid, while ¢, and y, are arbitrary solid harmonics of degree n. 
The stress exerted in the x-direction across a spherical surface 
of radius 7 by the fluid outside the sphere is p,., where 


Nl A0.0s 2n?+4n+3 oms 9 Dn | 


Pra = 2 be "Sa! Mtl) QnFl) Qnt8) de re 
(3) 


0 hn G) Xn 0 Xn 
m— 1 — — ——— ps 
+nd ¢ | ee 
The corresponding expressions for rp,, and rp, are obtained by cyc- 
lical interchange of x, y and z. The total force exerted on the sphere 
in the x direction is given by the integral over the sphere surface 


F= Sf Dreds. (4) 
It is evident from the orthogonality properties of surface harmonics 
that only those in (3) of zero order survive the integration, and fur- 
ther consideration shows that all of these drop out except the term 
E x/r*> in p-.. The final result is 


at Es (5) 


so that in any such fluid motion the total x-force on a sphere depends 
only on the coefficient of x/r* in the expansion of the pressure field 
in solid harmonics about the center of the sphere as origin. Similarly, 
the total y-force on the sphere depends only upon the coefficient of 
y/r*, and the z-force only upon that of the z/r? term. 
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II 
Such a fluid motion as considered above is uniquely determined 
throughout a region by its velocity components over the boundary 
(Lamb, p. 584). We proceed now to determine the value of E in terms 
of the boundary velocities on two concentric spherical surfaces of radii 
pand J . 
To this end we consider the quantity 1 defined by 


L=xu-+ yv + zw, (6) 


which is simply r times the outward radial velocity of the fluid. In 
terms of the above quantities its value (Lamb, p. 564) is 


i nr2 


= 7 = Danps) Pn FEM on - (7) 


: 2(2n+8) 


In this expression E' appears only in connection with the first order 
surface harmonic «/r. The following quantities are thus involved 
with # in the boundary conditions on 1 ; 


D,: Ax 
D.g2% Os Peo 
oi: Hx (8) 
g-2? Ga/r 
Upon equating coefficients of the x/r surface harmonic we obtain 
Ls eucnithesid Tied ii 
072 da: pe OG ae be 
(9) 
rd te oe 1 
ion vay adhered cy aoa 
where 
3 
= 10 
fi; Teste SS «lds (10) 


is 1/r times the coefficient of x/r in the expansion of / in surface har- 


- monics on a sphere of radius 7. It is the average over the sphere sur- 
face of 3 times the x component of the outward radial velocity of the 
fluid. 

Similarly, considering the boundary values of the velocity com- 
ponent w and picking out the terms involving zero order surface har- 
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monics gives 


A 2E 1 
ah © pol Riess 8 id 22 1-0 
ei T 35 t = 2 
(11) 
A Zia Oo 
cr Tg a Uy , 
where 
1 ; 
oS ee 12 
Ur Te uds (12) 


is the coefficient of the zero order surface harmonic in the expansion 
of u over a sphere of radius 7. It is the average over the sphere sur- 
face of the z-component of the fluid velocity. 

Solving from (9) and (11) 


OT (T5—t°) Au —5T (T?—1?) (Ayu—Al) 


5 Tt(P?—t)? — 4(T—_2) ("—P) (13) 


E=3yt 


where 
A = 1 


Age = Tp — OU; (14) 
AS Te 


Thus the x force on a sphere due to any such fluid motion is ex- 
pressed in terms of the quantities (10) and (12) which are certain 
averages of the fluid velocities over the surfaces of two concentric 
spheres. If v = w = 0 on these surfaces, while w assumes constant 
values, then 4,u = Al and the result as expressed by (5) and (13) 
reduces to that of Cunningham: if then T is made much larger than 
t it goes over into Stoke’s formula. 

This work was aided in part by a grant from the Dr. Wallace C. 
and Clara A. Abbott Memorial Fund of the University of Chicago. 
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PHYSICOMATHEMATICAL ASPECTS OF SOME PROBLEMS 
OF ORGANIC FORM 


N. RASHEVSKY 
THE UNIVERSITY OF CHICAGO 


: In connection with previous mathematical studies on cell polar- 
ity, the possible application of the results obtained before to differ- 
ent embryological phenomena is discussed. Methods for a quantita- 
tive mathematical approach to such phenomena as gastrulation, for- 
mation of different folds, closing of a half blastula, etc. are outlined. 


Of all the problems of biology, the problem of form and struc- 
ture of multicellular organisms seems to be the most difficult from the 
point of view of a physicomathematical approach. Even the highly 
complex phenomena observed in the functioning of the central ner- 
vous system and covering the fields of neurophysiology and psychol- 
ogy, have recently been brought within the scope of mathematical 
treatment (Rashevsky, 1938, 1940a; Landahl 1938, 1939 a, b, 1940; 
Householder 1938, 1939, 1940). In spite of the amazing progress of 
experimental embryology and a discovery of a very large number of 
important and entirely new facts, the complexity of this field seems 
to increase with increasing study. The quantitative element is still 
almost completely lacking in the experimental researches, and this 
seems to be due not so much to a lack of efforts on the part of the 
experimenter to introduce quantitative measurements, as to the elu- 
sive character of the phenomena studied, which makes it almost im- 
possible to point out where and how any quantitative method can be 
applied and what actually can be exactly measured. 

Under such conditions an attempt to present a completed and 
systematic theory of these phenomena would naturally be premature. 
But like in all branches of theoretical science, the theory of a group 
of phenomena must of necessity be preceded by a number of prelimi- 
nary, abstract theoretical studies, which deal with general conceiv- 
able explanations of the phenomena and study a large number of pos- 
sible explanations, although only one of these possible cases will per- 
haps be eventually applied (Rashevsky 1938 a, b; 1940 a). In this 
early stage any theory is of necessity incomplete and fragmentary, the 
fragments being rather disconnected. But even so, such a theory may 
at least suggest, where to look experimentally for quantitative rela- 
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tions, though it may not as yet give any consistent and adequate ex- 
planation of known phenomena. This paper is a contribution to such 
an early stage of the theory of organic form, and the reader is there- 
fore specially requested to keep in mind the above mentioned limi- 
tations. 


I 

In two previous papers (Rashevsky 1940 b, d), we have shown 
how under relatively simple conditions a cell may possess an auto- 
matically regulated axial polarity. This polarity is still of a very sim- 
ple type of symmetry, much too simple to be applied to any real cases. 
Nevertheless it possesses some interesting properties, one of them 
being that when such a cell is in any way divided in two parts, each 
part restores the same type of polarity, which is possessed by the 
original cell, regardless of the way in which the division is made. This 
latter statement however requires a qualification which, while limit- 
ing it to some extent, yet opens at the same time some interesting 
possibilities. 

The polarity, which consists in an axial gradient of concentration 
of a metabolite from pole to pole, may be conveniently measured by 
the difference «* of the average concentrations in the two hemispheres. 
This difference is given (Rashevsky 1940 b) by the approximate ex- 


pression: 
~..2 |8(Aaan—1) 
at = = | ae ia cach (1) 


NV Adan 


The symbols a, a and n have the same meaning as in the preced- 
ing paper. As regards the constant A, it has two different expres- 
sions, depending on the type of metabolic reaction considered. If we 


consider the case of a cell producing a substance, which diffuses out- 
ward, then (Rashevsky 1940 b): 


_ r2(2D + rh) 
A= 3D QD Brey ’ a 


the notations being again the same as before. If on the other hand we 
consider the case of a substance, which is produced in one part of the 


cell and consumed in another part, and does not leave the cell at all, 
then (Rashevsky, 1940 d): 


A= 
720 - 8D pe 


b being a reaction constant, defined in the previous paper. 
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The variation of A with radius 7 of the cell is quite different for 
the two cases. In the case represented by equation (2) A increases 
from zero to infinity, varying for large values of 7, like 7,2. In the 
second case, represented by equation (3), A tends to an assymptotic 
value 1/b. For such values of 7), for which 7,2b >> 3D, or 7,2 >> 
3D/b , the quantity A is practically constant. Taking b ~ 10° sec? 
and D ~ 107 cm? sec, we find that A will be approximately constant 
for r, > 102 cm. 

Let us first consider the case corresponding to equation (38), as- 
suming that 7 is within such a range as to make A approximately 
constant. Let us now divide the cell into halves equatorially. For 
- each half the quantity «* will again be determined by equation (1). 
But while the quantity A remains the same, has now changed, be- 
ing different in the two halves. For n denotes the averuge concentra- 
tion of the catalyst in the cell, and according to the discussion of the 
two preceding papers (Rashevsky 1940 b, d), that concentration is 
different in the two hemispheres. In the notations of the previous 
paper, the concentrations in the two hemispheres are nm, and n,, and 
the difference n. — 7, is given by equation (23) of loc. cit. (Rashev- 
sky 1940 b). On the other hand, x is given by equation (27) of the 
same paper. It follows therefore, since x* is the root of that equation 
(27), that 


* 


1 ah pee (4) 
Combining equation (1) and (4) we have: 
"3 (Aan —1) 
Nz — TM, = 2 \ ee metry. Ce oe (5) 
On the other hand we have (Rashevsky 1940 b): 
N+ Ny, = 2n. (6) 
Equations (5) and (6) give: 
3 (Aaan — 1) 
=a -_/ SS — Ataoen an (7) 
and 


r 3(Aaan — 1) 
age al ao CT ENE ey 


When we divide the original cell equatorially into halves then, 
while in each half x* is given by equation (1), yet the role of m in 
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one half is taken over by n, < 7 and in the other by m > n. Since, 
by the nature of its derivation, equation (1) holds only for Adan > 1, 
and at the same time when Aaan is only very slightly greater than 
one, an increase in ” results in an increase of «*. Hence, of the two 
half cells, the one that corresponds to pole N° 2 of the original cell 
will have a quantitatively larger polarity than the one corresponding 
to pole N° 1. 

As regards the half cell corresponding to pole N° 1, it can be 
shown, that with the particular mechanism of polarity adopted here, 
and leading to equation (1), it will loose its polarity after the division, 
for 7, will decrease to such an extent that Aaan, will become always 
less than 1. A somewhat different, more complicated mechanism 
which shall be described in a separate paper, can however be con- 
ceived, in which this limitation does not hold, and in which both 
halves preserve their polarity. In one half however the polarity is 
enhanced, in the other weakened. 

A further division of each half cell again in an equatorial plane 
will result in four quarter cells, in which the polarity will increase 
from one extreme cell to the other. Thus if an originally polar cell is 
divided by a series of equatorial sections into a group of smaller cells, 
the polarity in the cells of the group originating from pole N° 1 de- 
creases, that in the other cells — increases. But as n decreases more 
and more in one half of the cells, a situation will be reached when 
Aaan < 1. From now on the cells for which this situation obtains 
loose their polarity entirely, while the polarity of the other cells in- 
creases. If, as in the segmentation of an egg, the division into 
daughter-cells is not accompanied by any appreciable growth, the 
daughter cells will become smaller and smaller, until finally their 
radii 7; become so small that 7; < 8D/b . In that case A begins to de- 
crease. The decrease of A enhances the decrease of «* in the cells 
originating from hemisphere N° 1, and slows down the increase of 
x* in the cells, originating from the hemisphere N° 2. As 7; , and 
therefore A, tends to zero, Aaan —1 decreases more and more and 
thus «* will reach a maximum and then decrease. 

We now have the following situation. If an originally polar cell 
is divided into a few parts, each part preserves the same qualitative 
polarity, although quantitatively the parts will be different. As the 
segmentation proceeds, a process of differentiation takes place in the 
aggregate of cells. At one pole of the aggregate the cells gradually 
loose their polarity, while at the other pole the polarity of the cells 
increases, and if as a result of segmentation, the individual cell be- 
comes sufficiently small, it reaches a maximum and begins to decrease. 
An equatorial section of such an aggregate in two, at the proper stage, 
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results now in two qualitatively different halves: one polar, the other 
non-polar. 
While this picture is highly oversimplified, it offers nevertheless 
a model of a phenomena commonly found in embryology. Moreover it 
gives us a quantitative picture of what happens. 
Quantitatively the change of the average n for each division is 
given by the non-linear difference equations: 


3(A;.dan';.. — 1) 


N= Nia — 
As,07030' 5-4 


(9) 


for the pole N° 1, and 
5) ees Nee (10) 


3. _g302n". 
jaan ia 


for pole N° 2. For the case of equal division the value A; is defined by 


(0-8) ?? 7? ; 
(0-8)?* 67)? + 3D’ 


where 7, is the initial radius. The corresponding values of x* are then 
obtained from equation 1 into which the values of 7’; or 1’; and of 
A; are substituted. If the rate of segmentation is given as a function 
of time, then 7 is also a function of time, 1 = f(t), and the whole 
process of differentiation is described quantitatively in its dependence 
on time. A direct measurement of the polarity as defined here may 
be impossible. But by studying different theoretical cases of the in- 
fluence of polarity on other properties of the cell, as we shall do in the 
subsequent sections, we may derive relations between directly observ- 
able quantities. But even these considerations suggest the measure- 
ments of some quantities which may be used as a natural measure of 
polarity. 

A more complex situation may also be considered, by assuming 
that the quantity of catalyst in the system does not remain constant, 
but varies with respect to time according to some definite law, ex- 
pressed by a differential equation governing the kinetics of the reac- 
tion responsible for its production or destruction. 


Ai= (11) 


II 
In a previous paper (Rashevsky 1940 c) we have seen how un- 
der the influence of different forces, the cells of an aggregate, formed 
by the segmentation of a single egg, may arrange themselves in the 
shape of a hollow spherical shell. Combining these results with the 
considerations of the preceding section we find that in such a shell 


114 MATHEMATICAL BIOPHYSICS 


the cells in the neighborhood of one pole will be relatively strongly 
polarized, whereas in the other part of the shell they will not be Ppo- 
larized at all (Figure 1). This results in an asymmetry of the physical 


FIGURE 1 


properties of the shell, and it is of interest to study possible effects 
of such an asymmetry. To this end we must consider what effects a 
polarity of the nature discussed here may have on other properties of 
a cell. Of the many possibilities we shall consider here only three. 
a). The difference of physicochemical properties due to the dif- 
ferences of concentrations at the two poles may result, among other 
things in a differential growth of the cell at the two poles. Consider 
for simplicity a plane disc composed of one layer of cells, polarized 
in the direction normal to the disc (Figure 2). If due to the polariza- 


Uebedets | 


FIGURE 2 


tion the upper surface of the disc grows more rapidly, than the lower 
one, then the disc will bend, becoming convex upwards. Let the origi- 
nal diameter of the disc be equal to S. Let due to differential growth, 
the diameter of the lower surface become S, and that of the upper sur- 
face S, > S,. Let S. — S, = 6 and denote the thickness of the disc 
by A. Considering that the difference in growth rate is the same for 
all cells of the disc, we shall find that it will assume a shape of a 
spherical sector. Denoting by # the angle which is subtended by such 
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a bent disc from its center of curvature and by r the radius of curva- 
ture of the lower surface, we have (Figure 8): 


6= ¢A, (12) 
and 
h = r(1—cos$) (13) 
Furthermore we have 
= Sa ; (14) 


Substituting the value of ¢ from equation (12) into equation (13) 
and (14) and substituting the value of r from equation (14) into 
equation (13) we find: 


ees 


te ant 


(1 — cos us 2 (15) 


If we now have the rate of growth of both surfaces given as func- 
tions of the polarization «* and if the latter is itself prescribed as a 
function of time by consideration of the preceding section, then S, , d 
and 6 are known functions of time, and equation (15) gives us b, the 
“depth” of the buckling disc as a function of time. 

As an illustration consider the case that 


dS ZU dS. 


— * oe = Ko e 
dt Ax ; dt AgX ? (16) 
M—d,=—a>0, 
while A remains constant. 
Then, if 
z* = u(t) (17) 
we have 
. S,=Sta, ftu(t)dt =S+aU(t) ; 
S.=S+a, fiu(t)dt =S + aU (t) ; (18) 
oa 0U (tT). 
Hence, introducing equations (18) into equation (15), we have: 
[S + a,U(t)]A4 aU (t) 
= = : 19 
iu aU (t) eS a ya 


As can be readily seen by expanding cos (6/24), equation (15) gives 
h = 0 for 6 = 0. The value of h reaches its maximum in the neigh- 
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borhood of 6/24 = 2, which corresponds to ¢ = 2a, ora completely 
closed shell. For 6 = 0 we have h < 0, which means a disc convex 
downward. 

If due to the polarization of the cells in the shaded region of Fig- 
ure 1 the inner side of the cells grows faster than the outer side, that 
region will buckle inward, and we shall have a phenomenon, remind- 
ing of gastrulation. Equation (19) gives us the depth of the invagi- 
nation after it actually started, as a function of time. Making differ- 
ent assumptions about the effects of the polarization upon the differ- 
ential growth rate, we shall have different forms of equation (19), 
which may be compared to different observed cases. 

True enough, the invagination in a gastrula does not proceed ac- 
tually in such a simple way. But these simple considerations should 
pave the way for more complex cases. By considering different dis- 
turbing factors, such as nonuniform thickness of the blastula, etc., 
we may arrive at a mathematical expression of more complex forms 
of invagination, which come closer to the actually observed ones. It 
must be remarked however, that while for instance in Amphibia the 
process of invagination is very different from the simple scheme here 
considered, in some cases, like Arbacia or Amphioxus, it is, while 
still complex, to some extent comparable with our oversimplified pic- 
ture. In any case the importance of having the phenomenon of invagi- 
nation properly timed, describing the depth of invagination quantita- 
tively as a function of time, becomes apparent. Hitherto no observa- 
tions of that kind have been made. In order to obtain smooth results, 
it will probably be necessary to use a very large number of observa- 
tions, and take averages, due to individual variations. But in prin- 
ciple here is a definite quantitative problem of embryology, which 
awaits a quantitatively minded embryologist to tackle it experimen- 
tally. 

That actual gastrulation is probably due to some such process as 
above discussed, seems to be made rather probable both by the model 
experiments of J. Speck (1918) as well as by the recent observations 
on the invagination of the isolated vegetative parts of Dendraster 
excentricus by A. R. Moore and Agnes S. Burt (1939). Ina previ- 
ous paper (Rashevsky 1940 e), we have considered asymmetric dis- 
tributions of metabolism in a blastula. If the effective rate of meta- 
bolism (Rashevsky 1940a) is positive everywhere in the blastula, ex- 
cept in a region near the vegetative pole, then the diffusion forces in 
that region will be directed inward, and cause that region to invagi- 
nate. Taking these forces to be of the usual order of magnitude for 
diffusion forces (Rashevsky 1940 a), considering that the principal 
resistance offered to invagination is the flow of the blastocoel liquid 
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through the invaginating part, as the blastocoel gradually closes, and 
taking for this resistance a value used previously (Rashevsky, 1938, 
p. 106), we find for the time necessary for a complete invagination a 
value of the order of magnitude of 10 hours, which agrees with some 
observations (Taylor 1937). From this picture one could also derive 
a relation between depth of invagination and time. Such a mechanism 
is however made rather unlikely by both the above mentioned experi- 
ments of Moore and Burt and by those of J. Brachet (1934). In the 
former an invagination is observed in open isolated parts of blastula, 
in which the conditions of a closed metabolizing shell are definitely 
not satisfied. In the latters, it has been shown that gastrulation may 
proceed normally in the absence of oxygen. Since oxygen is the only 
substance with negative rate of metabolism, used by an early embryo 
in quantities sufficient to produce appreciable forces, Brachets’ ob- 
servations seem to rule out the above possibility. Thus some sort of 
polarization resulting in an autonomous invagination of the vegeta- 
tive region, regardless of whether it forms a part of a whole blastula 
or not, seems to be the cause of gastrulation. 

b). We now consider a different possible mode of action of cell 
polarity. 

Referring again to Figure 2 we may consider the case, that due 
to a difference of concentrations on both sides of the disc, the surface 
tension on the upper side is smaller than on the lower. This will 
again result in a convexity upwards. The calculation of the variation 
of the convexity, with time would proceed as follows. We first calcu- 
late the buckling of an elastic disc under the influence of such forces. 
Then apply the analogy between elastic deformation and plastic flow 
(Rashevsky, 1938 a, 1940 a). 

c). The most interesting consideration however is perhaps one 
which does not involve any additional assumptions. Considering again 
Figure 2 and assuming that the polarization is of the type correspond- 
ing to equation (3) and is so directed that the concentration is higher 
at the top surface, than at the bottom, we have a situation represent- 
ed on Figure 4. If c: > C2 > ¢s, then, since c plays the role of a po- 


FIGURE 4 
tential, the surface forces are distributed as indicated by the arrows. 
This again results in a buckling upward, the rate of which may be 


calculated in a similar manner as suggested above (C.f. also section 
III). In this way we would have an equation for the rate of invagina- 
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tion, expressing the latter in terms of essentially the same constants, 
which have been used previously in cell biophysics. 


Ill 

If the origin of cell polarity may be considered as being connected 
with the axial gradients of concentration, then a spherically sym- 
metric cell may be said to have radial polarity. The axial polarity, as 
studied here is due to the interaction of a metabolite with some cata- 
lyst. By considering two or more independent metabolites of which 
one behaves in such a way as to produce an axial polarity, while the 
other has an approximately spherical distribution, we obtain a cell in 
which a radial polarity is superimposed upon an axial one (Figure 5). 


FIGURE 5 


In fact such a superposition may occur even in case of a single meta- 
bolite. In the case studied before (Rashevsky 1940 b) which corre- 
sponds to equation 2, the concentration of the metabolite is smaller at 
the periphery than near the center, only the maximum of concentra- 
tion lies excentrically, giving thus an axial polarity. 

As a result of the segmentation of such a cell, we shall have a 
blastula, which also exhibits the mixed polarity. If we consider the 
radial component and its effects, we see that by a similar argument as 
used before, we may make the radial component responsible for the 
closing of a half blastula, obtained by cutting a whole blastula in two. 
By similar reasonings as before we may obtain an equation giving us 
the rate of closure and its dependence on the initial conditions, that 
is on whether we observe a half blastula or a smaller or larger frac- 
tion of it. The phenomenon of closure of parts of blastula is well 
known (Horstadius, 1939; Weiss 1939). However only qualitative 
descriptions of the phenomenon are available. Quantitative observa- 
tions, which in this case would not be so difficult at all, would give us 
valuable information as to the mechanism responsible for the closure, 
for different mechanisms will lead to different relations for the de- 
crease of the opening with respect to time. 

Compared with the problem of invagination, the problem of clos- 


N. RASHEVSKY 119 


ing of a part blastula presents an additional interesting feature. Con- 
sider a partly closed blastula (Figure 6), whose cells do not possess 


Ca 


Ci 


C3 
FIGURE 6 


any intrinsic polarity and into which a substance diffuses from out- 
side. As can be readily seen from general physical considerations and 
proven by the application of the usual approximation method, when 
the external diffusion coefficient D. is finite the concentration c, at 
the inner surface will be smaller than the concentration c, at the outer 
surface, so that c, < c.. Thus inside the cellular layer the cells will 
now show a polarity, but this polarity is so to say an induced one, 
persisting only as long as the cell is a part of the whole aggregate, 
having the shape of a non completely closed shell. The polarity will 
be the smaller, the more “open” the blastula and will vanish complete- 
ly for a perfectly flat disc. Considering now the effects of this polar- 
ity either on differential growth or on surface tension or directly upon 
the diffusion forces, we shall again obtain different expressions for 
the rate of closure, which may be compared with experiments. 

The above discussed “induced” polarity raises another important 
question, namely: if an intrinsically polarized cell is brought into an 
external diffusion field, will this field tend to orient the polarity of 
the cell in a definite direction? If yes, will the forces exerted by the 
external field upon that cell be such as to orient the whole cell, or such 
as to change the orientation of polarity within the cell? This prob- 
lem, which may present some interesting physicomathematical as- 
pects, is of interest in connection with the experiments of implanting 
additional micromers into Arbacia blastulae. 


IV 

We shall now consider still another possible effect of polarity. 
Consider now a closed shell of thickness 4, consisting of several lay- 
ers of cells (Figure 7). Such a multicellular layer possesses a certain 
degree of elasticity like any other tissue. Let those cells be polarized 
radially. Denoting by 71 and r, the inner and the outer radii, by S; 
and S, the length of the great circles on the internal and external sur- 
faces and putting S. — S: = 6, we have 
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FIGURE 7 


Se— Si 6 
= — SSS 2 
A = Foe ty = a’ (20) 
If 5, due to polarity, increases with time, so that 6 =U (t), then the 
thickness A also increases. This results in radial tensions inside of 
the layer, the average radial stress being given by 


A — Ay 
Ae ea 


where E denotes the elasticity modulus of the tissue, and A, the in- 
itial thickness of the shell. We may speak approximately of the tensil 
strength o; of such a tissue. When a, > o; , the tissue will break in the 
radial direction and form approximately concentric layers with a 
space in between. If 6 = U(t), then A = U(t) /2a, and the time at 
which the rupture, or speaking in embryological terms, delamination 
will occur, is given by 


on, — E (21) 


U(t) — 2a Ay 


ae OF ewes (22) 


If an axial polarity is superimposed upon the radial one, so that the 
radial polarity varies from point to point, being say, largest at the 
point P, (Figure 7), a smallest at P., then 6 is a function of the posi- 
tion on the surface and may be conveniently described as a function 
of the angle ¢ (Figure 7). Thus. 


a= Us, t), (23) 
and equation (22) now becomes 


U(g,t) —2QaAy _ 


. aA, ys (24) 


Solved with respect to t, it gives t as a function of d, 


t=v(¢). (25) 
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Thus now the layer splits into two lamina, gradually beginning with 
one pole and ending at the other. While the simple picture used here 
does not correspond to any real situation, yet it does remind of the 
delamination of mesoderm in Amphibia, and definitely suggests a 
quantitative study of the process of delamination with respect to time. 

The author is indebted to Mr. H. D. Landahl for checking the 
calculations. 
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The average concentrations of a substance diffusing into or from 
an open spherical shell, in which it is consumed or produced at a 
constant rate, are calculated by the approximation method. An ap- 
plication of the result to the problem of deformation of such a shell 
under the influence of diffusion forces is indicated. 


In the preceding paper (Rashevsky 1940) we have indicated the 
interest of studying the diffusion processes and mechanical forces in 
cellular aggregates, which have the shape of open spherical shells. 
In the present paper we shall give as an illustration a very crude 
treatment of the problem of deformation of such an aggregate. 

Consider an open spherical shell, the meridional section of which 
is represented on Figure 1. Let a substance be metabolized inside the 
body of the shell at a constant rate q gm- cm“ sec. For simplicity we 
shall consider the permeability h of the surfaces of the shell as infinite. 
We shall denote by D; and D, the internal and external diffusion co- 
efficients respectively. The average concentration inside the shell we 
shall denote by é, while ¢,, c. and c, denote the average concentra- 
tions at the different surfaces of the shell, as is clear from Figure 1. 
By c. we shall denote the average concentration in the plane of the 
“opening” or “mouth” of the shell. The remaining notations can be 
understood from Figure 1. We shall put 7: + 6=7,. We assume 
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6 << 7r,, and shall therefore put everywhere 7, = 72 = 7, except 
when 7, and 72 occur as a difference 72 — 71. 

Denoting by S; and S. the areas of the internal and external sur- 
faces respectively and by S; the area of the surface formed by AB ; 
denoting further by V the volume of the shell, we have 


S, = 2ar,2(1— cos ¢) ; S2=2ar.?(1— cos ¢) ; 
S,=2ar'6=2aridsing; ”=7,sin¢g ; (1) 
S,=ar2=ar,? sin’?¢ ; V = 2ar?(1— cos ¢$)6. 


The concentration c inside the shell has its average value € somewhere 
between the dotted line and each surface, so that the average gradi- 
ents inward and outward are of the order of 4(é — ¢,)/6 and 4(é — 
c.)/d. From physical considerations it can be seen, that the gradient 
near the opening and in the direction of increasing angle ¢ is of the 
order 4(¢ — ¢c;)/6. Considering as before that outside of the shell 
the gradient extends over a region of the order of 7 , we now have the 
following equations for the material balance: 


€ — Cz 
6 


Moreover we have a flow inside the shell, due to the fact, that c, is 


different from c,. The average distance 6 over which the difference 
c, — c, extends is of the order of + O’ O” (Figure 1), hence 


€— Cz Cz — Co C3 — 6p 


=D. > AD; 
fh: 


4D; =D, (2) 


5 = 5 (1—cos¢) . (3) 
Hence 
C—C, C1 — Cy 
$Daraayer tae Dineey at taarayt a 


Moreover the total flow through the surface S, equals the total flow 
through the surface S,. Hence 


2D.S, Cy Cz C4 — Co : 


r(i—cosgy 


(5) 
Finally we have the equation expressing the equality of the amount 
metabolized with that flowing out: 


€— Ce 
Fy 


€— 


ré6(1—cos ¢) q = 4D; 5 


r(1—cos ¢) + 4D; 


(6) 
Xx r(1—cos ¢) + 4D; (é—e¢,) sing. 
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From equations (2), (4), (5) and (6) we can determine the quan- 
tities €, €,, C., C3, C, by elementary, though somewhat cumbersome 
calculations. Putting 


sin¢g=—s,; 1—cos¢=&, (7) 
we find 
€=Co+ 
(8) 
OS82qQr (4Dir + D.6) (8Dirs, + 2Di1rs,*s2 + D.68;") 


4D;D-. (8D;7s_ > 2D 7S77S85 = D-.6s?;) (7s; + 68,) + (4D;7r + D.6) 781782 


When g>0, then €—c,>0, as should be the case. For d6=a, 81 
= 0, s, = 2, the shell becomes closed. Then 


ERM yea g prem PM (9) 


an equation, that can be obtained directly. 
We shall be interested in the quantity c. —¢, ; this is given by 


2D;D.76 (2812 — 482 — $1782) (€ — Co) 
C — & = — OO TOOT ae (10) 
(4D;r + D.6) (8Dirs, + 2Di781?S2 + D-68,") 


It can be shown that the expression (28,2 — 48, — $1782) is always 
negative. Therefore when €—¢,<0, ¢:— > 0, or when g<0, 
then c.—¢i:>0. 

For q < 0 the shell is subject to an inwardly directed pressure 
of the order of magnitude of 


p =" T(@—a) - (11) 

which tends to close the shell. To find the equation governing the 
process of such a closure, we consider the corresponding elastic prob- 
em. An elastic shell of thickness 6, subject to a uniform pressure p is 
deformed so that the radius 7’ of the opening changes by an amount 
given by (Hetenyi 1939): 


Ar’ pr 

— = — py); 12 

where E is the elasticity coefficient, and u the Poisson’s ratio. 
Consider the shel! composed of cells as being a plastic body of an 

average viscosity 7 (which is not to be confused with the viscosity of 

protoplasm of the cells composing the shell), we find (Young 1939) 


the differential equation for the variation of r’ with respect to time by 


Se 
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putting in equation (12) w=1/2, E=7, and substituting dr’/dt 
for Ar’ . This gives: 
Wipe eS em pr (13) 


By means of the fourth and sixth of the relations (1), we may ex- 
press 7’ and r( =—7,=—7;) in terms of ¢ and V . Considering 6 as con- 
stant during the deformation, we thus have a differential equation for 
the variation of ¢. Such an equation could be compared with the rate 
of closure of a half blastula, if quantitative data were available. 

The author is indebted to H. D. Landahl for checking the calcu- 


lations. 
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_ A general solution of the formal nerve conduction problem is 
given. As illustrations of the general method, the capacitative sin- 
gle-factor and the non-capacitative Lapicque problems are solved. 
Comparisons between velocity formulae for capacitative and non- 
capacitative models indicate that previously determined non-capaci- 
tative velocities are considerably too high. 


The formal theories of nerve conduction as presented by N. Ra- 
shevsky (1938, part Il) and W. A. H. Rushton (1937) have usually 
included the simplifying assumption that the network representing 
the nerve axon is non-capacitative. An attempt to extend one of Ra- 
shevsky’s theories to include capacitative effects has been made (Wein- 
berg, 1939) but was rather unsuccessful, primarily because the prob- 
lem was treated from the integral instead of differential standpoint. 

In the present paper a general solution will be outlined for all 
Rashevsky-type formal nerve conduction problems; in principle, by 
the use of this method the steady state velocity of conduction can be 
calculated for any generalized capacitative nerve network in which 
excitation is governed by any one of the formal excitation theories 
(Rashevsky, 1938, part II). As an illustration of this method, the 
velocity of conduction will be calculated for some simplified models 
among which may be found all those considered by previous writers. 

General Equations. The formal nerve conduction problem is di- 
vided into two separate parts: first, on the basis of some assumed 
electrical structure the distribution of currents in the nerve is calcu- 
lated; and second, the currents so calculated are introduced into an 
excitation equation which describes the time variation of excitatory 
state as a function of the stimulating current. From the condition 
that excitation occurs when the excitatory factor exceeds a threshold, 
the velocity of propagation may be calculated. 

Suppose that some general electrical network has been assumed 
to represent the nerve, and suppose the steady state has been reached. 
Then the differential equation satisfied by the local bio-electric cur- 
rent, I, will be of the form 


F(a" 1/dé*, dv I/dé"* ,--- 1) =0, (1) 
127 
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where € = x—vt measures the distance x from an origin moving 
with the disturbance whose velocity is v; equation (1) is obtained 
from Kirchhoft’s laws in the usual manner (Webster, p. 44). The 
solution of equation (1) is 

(ees (2) 
whether this solution can actually be calculated in closed form de- 
pends, of course, on F — that is, on the the electrical character of the 
assumed network. 


The excitation equation governing the variation of an excitation 
factor y, will in general be a differential equation of the form 


f(d*y/din ,d™1y/dt™...y, dl /dt, d-1/dt™ -.. 1) =0 (3) 


which may be solved for I as a function of y and the derivatives of 
y and I: 


I= wy, dy/dt ,--- dI/dt, a? /dt? ---), (4) 


or, in the steady state, since d/dt = —v(d/dé) , 
1 = p(y,—vdy/dé, --- —vdI/d&, v2d?l/d& «.-). (5) 


The derivatives of I can be calculated from equation (2), and if they 
are substituted into (5), y becomes a function of y, &, and v alone. 
Upon equating the two expressions for J given in equations (2) and 
(5), one obtains 


o(€) =yply, —vdy/dé , v?d?y/d&---), (6) 


which is a differential equation in y with solution y = y(é,v). The 
steady state velocity is then determined from the condition that at 
E=0 (ie., x = vt), y must just reach the threshold value y : 


Yo=y (0, v). (7) 


Equation (7) is an algebraic or transcendental equation in v which has 
the solution 


Vv = V(Yo) (8) 


and this velocity formula may be expressed in terms of the rheobase 
since Yo is related to the rheobase in a known manner. The method 
just outlined is similar in character to that used by Rushton (1937) 
and Umrath (1928) ; it is much simpler than the functional equation 
solution of Rashevsky. 

Applications: As an illustration of the preceding general con- 
siderations, the velocity formula (or rather the algebraic equation 
(7) determining the velocity) for a capacitative nerve model obeying 
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a general two factor excitation law will be derived. A very simple 
capacitative model, shown in the accompanying figure, will be assumed 
to represent the nerve. The capacitances C may be identified with the 
capacitance of the nerve membrane, p with the transverse leak resist- 
ance of the membrane, 7; and r, with the interior and exterior longi- 
tudinal resistances, and E with the action voltage, all line constants 
being measured per unit length. In the present formal model, F is 
assumed to be zero until action occurs, at which time FE suddenly 
attains a steady value E,.* 


FIGURE 1 


Before writing down the differential equation of this circuit, it 
will be well to consider which of the network currents — 1%, the cur- 
rent through the transverse resistance p, 7, the current into the con- 
denser, or i,, the longitudinal current — should be identified with the 
exciting current of equation (3). Fortunately, this question does not 
present any difficulty, for in the steady state, as may be easily verified, 
all three currents are proportional. For this reason, any linear com- 
bination of these currents may be used to represent the exciting cur- 
rent, the only effect on the excitation equation being to change the 
numerical value of the constants in it, without changing the form of 
the function f of equation (3). Inasmuch as the numerical values of 
these constants are obtained only from experiments in which the 
exciting current is the external stimulating current (Blair, 1932), 
the values of the constants to be assigned for re-excitation by action 
currents are not directly determinable, and these values are not sig- 
nificant for the present. 


* The general argument given here remains essentially unaltered if a non- 
explosive change in H is assumed to occur, so that in the steady state E would be 
some continuous function of £, say E(é). This function would be closely related 
to the observed action potential spike. The differential equation for steady state 
conduction would now be continuous throughout the é range, but it could not be 
integrated, in general, except by numerical methods. This type of non-explosive 
change has, been discussed elsewhere (Offner, Weinberg, and Young 1940). The 
assumption that E remains at E, even after the impulse has passed does not 


affect the velocity appreciably, as is shown in the previously mentioned paper. 
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If i, is identified with the local exciting current, then the differen- 
tial equation for ip may easily be shown to be 
od? i/au? — MpC di/at — Mi + 0H /oau? = (9) 
where M = r,-+r., and the subscript p has been dropped from %» : 
In the steady state, since 0/ot = —vd/0é , 0/dx = 0/0&, this equation 
becomes 


~ pd?i/dé? + vMpC di/d§ — Mi4+ @h/de = 0 (10) 
the solution of which is (Rushton, 1937 p. 216) 
Cages To alia (11) 


where i, is the action current through p at the junction of excited and 
unexcited regions, 1/L’ = [(MCv/2)?+ 1/L?]'-+ MCv/2, and L= 
Vp/M , the characteristic length. 

The most general two-factor excitation theory may be formulated 
as a second order differential equation 


dey /dt? + ady/dt + py = adi/dt + bi (12) 


where y is the difference between excitatory and inhibitory (‘‘accom- 
modation’’) factors and so measures the state of excitation of the 
nerve. Equation (12) is simply one way of writing the two first 
order equations of the Hill-Rashevsky excitation theory (Hill, 1936; 
Rashevsky, 1938, chap. 21); if b =0, the equation reduces to the 
“normal accommodation” case of the Hill-Rashevsky theory which is 
equivalent to the Monnier theory (1934). The constants a, 8, a and 
b are easily expressed in terms of the constants of the original 
theories ; however, for the present, it is sufficient to note that a and a 
are related to the rheobase, R, by R = ay,/a, where yp is the excita- 
tion threshold. 
Equation (12) may be written in the form 


Baebes: f ‘(dy /dt? + ody/dt + By — bi) at (13) 


as in equation (4) ; equating the two expressions for 7 from (11) and 
(13) and differentiating one obtains (for the steady state where 
d/dt = —vd/dé) 


dy /de* — (a/v) dy/dé + By/v? = (a/vL' + b/v?) i) es", = (14) 
The solution of this equation is readily seen to be 


_ (aL'v + dL”) ines!’ 
~  w Falot BL? * ite 
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Since at the junction between the two regions £=0 and y = Y, the 
algebraic equation defining the velocity [equation (7) of the previous 
section] is found from (15): 

big 

Yo 

Equation (16) involves high powers of v (since L’ is a function of 
v), and so it cannot be solved exactly, although numerical solutions are 
feasible. Certain simplifying assumptions can be made which will lead 
to equations which are more easily handled, and these will be con- 
sidered now. 

Capacitative Networks. The simplifications which can be intro- 
duced into the theory are of two sorts: on the one hand, the electri- 
cal network can be simplified by, say, eliminating the capacitative 
elements, and, on the other hand, the general excitation equations can 
be simplified in various ways. If the capacitative network is used, 
(16) is the fundamental equation defining the velocity for the general 
two-factor excitation theory. The various simplifications introduced 
into (16) by using simpler excitation theories may be listed as fol- 
lows: 

1. “Normal” accommodation: b= 0. 

2. Lapicque theory (Lapicque, 1926 p. 184; Rashevsky, 1931): 

B= O, 0520. 

8.. Single-factor, Blair (1932) theory: 6=6b=0. 

Case 8 is the one studied previously by the writer (Weinberg, 
1939) ; the asymptotic velocity for this case can now be calculated. 
Upon setting 6 = b = 0 in (16) and using the expression for the 
rheobase R = ay,/a, one finds the velocity to be 


vy=aL' (17) 


where 4 = (i, — R)/R is the “safety factor.” Equation (17) may be 
solved for v: 


an ( __ Wo 
tag 
L Yo 


Vv 
is ace ae =0. (16) 


ab alL 
V1-+ alMCL? © 


If, as Rushton (1937) assumes, a = 1/Cp, then equation (18) may 
be written as 


(18) 


sso Sears sean (18A) 
CVMp(1 +4) 


Non-Capacitative Networks. To illustrate further the power and 
generality of the present method, the results for a non-capacitative 
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network will also be considered briefly. For a non-capacitative net- 
work L'=L, so that from (16) the general velocity equation 1s 
(De ) bio 


es ay 7 tee ak 0 19 
au +b = 0, (19) 


which has the two positive solutions * 


adL 
2 


v= ( ) + V (aLA/2)? — L? (B — bt/Yo) . (20) 

The formula (20) includes all previous non-capacitative velocity 
formulas as special cases. For example, the Lapicque theory velocity 
which was shown to be a constant but whose exact value could not be 
determined (Rashevsky, 1931) is 


v = (alL/2) + V (ad /2)? + L*b%0/Yo (21) 


while the Blair theory velocity, which is essentially the same as the 
normally accommodated two-factor velocity (Rashevsky, 1938, p. 
212), 18 


v=aLi. (22) 


Since the Blair velocity is zero if the negative sign is used in (20) 
only the larger positive root is allowed there. 

A comparison of (18) and (22) shows immediately that if 2 can 
be assumed to be a physical quantity independent of the circuit con- 
sidered, the non-capacitative velocity, Unon-cap, iS greater than the ca- 
pacitative velocity, v.4,». An estimate of the difference between them 
can be made by assuming a © 1/Cp. Then from (18) and (22) 


Vnon-cap — V1 a A Veap 


and if A © 2, Vnon-cap © 1.7V cap , SO that the non-capacitative velocity is 
considerably larger than the capacitative. 

As has been mentioned already, the method outlined here is sufti- 
ciently general to handle any formal conduction theory in a straight- 
forward, simple manner. It seems unlikely however, that formal 
theories of nerve conduction will ever be able to explain the more 
intimate characteristics of propagation phenomena; rather, it would 
appear to the author that further theoretical development ought to be 
in the line of increasing the physical content of the theory, and, in 
particular, deducing the formal laws of excitation from structurally 
plausible, physical postulates. 


* As will be seen presently, only the larger root is of interest. 
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A NOTE ON THE HOROPTER 


ALSTON S. HOUSEHOLDER 
THE UNVERSITY OF CHICAGO 


‘ By assuming the fixity (but not the symmetry) of correspond- 
ing points on the two retinae, it is possible to derive the equation of 
any horopter when one is known. In particular when, as experiment 
shows, one horopter is linear, then all horopters must be conics. 
These have the form given by Ogle, but whereas Ogle leaves one 
parameter undetermined at each fixation, on our assumption the only 
arbitrary parameter is determined by the position of the linear 
horopter. 


It is well known, and easily verified, that when any fixed object, 
not too far away, is fixated by the two eyes, then near objects as well 
as distant objects appear double, whereas objects in some region 
which includes the object fixated do not appear double. That is, when 
the eyes are held fixed, the two retinal images of any given object will 
or will not be fused, depending upon its position relative to the ob- 
server and to the object he is fixating. Therefore it may be supposed 
that to each point on one retina there corresponds a unique point on 
the other retina, called its “corresponding point,” which is such that 
if the same object is imaged on these two corresponding points the 
two images will fuse. Thus there is established a one-to-one corres- 
pondence between the points on the two retinae, and if the optical 
system of the eye were specified, and the correspondence were known, 
then, the eyes remaining fixed, one could determine the locus of points 
in space whose images would fuse, i.e., whose images would fall on 
corresponding retinal points. This locus would be a piece of an ordi- 
nary surface in space, known as the horopter. 

Actually it is known that to a given point on one retina there 
corresponds a region of points on the other retina in the sense that 
fusion can occur, although the fusion becomes more and more difficult 
as the region is departed from. Certainly the boundaries of this region 
are not clearly defined. We may therefore think of each point on one 
retina as having a unique, “truly” corresponding, point on the other 
at which fusion is easiest, and that as this point is departed from in 
any direction the fusion becomes increasingly more difficult and the 
splitting of the images increasingly easy. We shall therefore continue 
to speak of retinal correspondence as a one-to-one correspondence. 
Images of the same object which fall upon points that do not corres- 
pond in this sense will (as is customary) be said to be disparate, even 
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though they may fuse for the observer. Also we may distinguish 
among nasal or medial disparity, temporal or lateral disparity and 
vertical disparity, the adjectives being self-explanatory. 

Depth perception may be “explained” in terms of the horopter 
by saying that objects whose images are medially disparate appear 
farther, those whose images are laterally disparate appear closer, than 
points on the horopter, this being true whether or not fusion has 
occurred. This, however, is not a true explanation, and for several 
reasons. First, it says nothing about the subjective form of the horop- 
ter and how it is achieved. Second, while we may grant that medially 
disparate images should be qualitatively different from laterally dis- 
parate images, whether or not these are fused, nothing is said about 
how these qualitative differences become organized into a unified, 
three-dimensional space with quantitative attributes (Roelofs, 1935). 
Finally, the manner in which fusion occurs remains a mystery. Never- 
theless, the notion of the horopter is a useful descriptive concept in 
discussions of depth perception. 

Let us restrict ourselves in the present discussion to the con- 
sideration of the horopter curve, which is the intersection of the 
horopter surface with the plane of the visual axes. We need to con- 
sider then only those light rays which lie in this plane, and this 
restriction will be tacitly assumed hereafter. Now it is convenient to 
suppose, first, that the light rays from any point to its retinal image 
can be represented by a straight line passing through some point in 
the bulbus — say the nodal point — which is common to all rays, and 
second, that the position of this point does not vary sensibly with 
respect to the retina as a result of accommodation. This is equivalent 
to saying that this common point can be regarded as the center of a 
pencil of lines through the points of the retina, that each of these 
lines gives the projection of a point in space upon a retinal point 
which is the optical image of the space point, and that the Euclidean 
properties of this pencil as determined by the retinal points are in- 
variant under changes in accommodation, as well as position. Where 
the image is not in focus on the retina we must regard the line as 
passing through the center of the circle of confusion on the retina. 

We now make a second assumption to the effect that the relation 
of correspondence, as defined above, is invariant under changes in 
position and accommodation. Thus if the retinal point P of the left 
eye corresponds to the retinal point Q of the right eye when the eyes 
are in one position, they will correspond when the eyes are in any 
position whatsoever. 

One attempt to give a rational derivation and mathematical rep- 
resentation of the horopter curve has resulted in the circles of Vieth 
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and Miiller, and is made on the basis of the further assumption that 
the eyes are essentially symmetrical, with corresponding rays (lines 
leading to corresponding points on the two retinae) making equal 
angles with the visual axes. These circles are therefore coaxal circles 
— the centers of the pencils being regarded as fixed in space — the 
common intersections of the circles being the centers of the pencils. 
Unfortunately, these circles fail to represent adequately the experi- 
mental results, for it is found that in all cases the empirical horopter 
has less curvature toward the observer, and in fact for distant vision 
the horopter actually has its concavity directed away from the observer. 
Lying intermediate to the horopters with concavity directed toward, 
and those with concavity directed away from the observer, there is 
found to be one which is a straight line (Southall, 1987; Ogle, 1938). 

The simplest family of curves possessing the properties here des- 
cribed is evidently a family of conics, and Ogle (1938) has derived 
the equations of the horopters on the assumption that they are conics 
passing through the centers of the visual pencils. For horopters 
which are symmetric with respect to the frontal plane these conics 
depend upon a single parameter, aside from those which are directly 
accessible, viz., the distance between the eyes and the location of the 
fixation point. This parameter, H, would be zero for the Vieth-Miiller 
circles, but empirically it is in general positive. Ogle’s derivation does 
not relate the values of H for different fixations of the same subject. 

This note is intended to point out that on the assumptions as 
to the invariance of the form of the pencils to the retinae and of the 
relation of correspondence, it is possible, given the horopter curve 
for one fixation, to predict the form of the curve for any fixation. In 
particular, from the fact that one horopter curve is linear, it follows 
easily that the other curves must be conics of the type described. This 
is true whether or not the linear horopter is perpendicular to the 
frontal plane. 

Let the centers of the pencils at the two eyes be located at 
(La, 0); let a and 2 — a be the inclinations (in the sense of elemen- 
tary analytic geometry) of the visual axes for the particular sym- 
metric fixation which yields the linear horopter ; let 6 be the inclina- 
tion of the linear horopter; let 6, and 6, be the angles made with the 
visual axes by a pair of corresponding rays (lines of the pencil passing 
through corresponding points on the two retinae). The assumption of 
invariance means simply that the functional relation between 6, and 62 
is independent of the position of the two eyes. Hence, since in this 
position the two pencils are in perspective (in the sense of projective 
geometry), the pencils remain projectively related, and therefore the 
the locus of their intersections is always a conic. 
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Any ray of the one pencil can be written in the form 
—extana+ (1+4/,)y—atana=0, 
and any ray of the other in the form 
xtana+ (1+1.)y—atana=0. 
The condition that the rays be corresponding is that they intersect on 
the line . 
xtanfB—y+atana=0, 


and hence that the determinant of the three linear equations shall 
vanish. This condition is equivalent to the condition 

4,/sin (a — 6) = 4./sin(a + 8) =», 
where » is an auxiliary parameter. Hence we have the equations for 
corresponding lines in the form 


—xtana+ [1+ sin (a—)ly—-atana=0, 


(1) 
xtana-+ [1—ysin (a+ £)]y—atana=0. 
Therefore 
tinepiee — sin (a — §) tan o 
sec?a + » sin (a — fp) ’ 
. (2) 
taiiia cee sin(a-+ #) tana 


sec? a—ysin (a+ 6) ° 
From this it follows that 
sin (a-+ £) cot 0, —sin (a —f) cot 0, =2cosacosB, (3) 
which is analogous to Ogle’s relation (10). In the special case that 
8 = 0, this relation becomes 
cot 3, — cot 3), = 2 cota. (4) 
On the basis of our assumptions, Ogle’s parameter H for the symmetric 
horopter is equal to 2 cota. 
Suppose, now, that the left eye is rotated through an angle @,, 


and the right eye through an angle w.. For symmetric fixation these 


angles are equal and opposite. The equation of the corresponding lines 
are then 


tan (a+ Shp fan ort sin {ore Ban eased 


1 — tan a tan w, + y sin (a — p) cta’ 
and 


Ne ee See 
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ade peeney ee nea le —tan a+tan w.—-» sin (a8) tan @2 _ Y 
1+tan a tan w,—yv sin(a+f) x—a,’” 
or 


(z+ a) (tana tan w,) — (1— tana tan w,)y +» sin (a — f) 
Xi(%#-+ a)tan o.— y] = 0, 
(5) 
(*« — a) (tan ao — tan wm) + (1 + tanatan w.)y + »sin (a + f) 
y PaO. 


By eliminating » between these two equations we obtain the equation 
of the locus of intersections of corresponding lines which is the horop- 
ter. This can be written in the form 


X [(#@—a) tan ow. 


(a +a)sin(a+@,)—ycos(at+@) [(«#+a)sinw,—ycosm,]sin (a—f) 


fia)... | (4-1) Sina. 
=0. (6) 


Thus when the linear horopter is known, the horopter corresponding 
to any other fixation point can be determined. 

The horopter is symmetric when and only when the coefficient of 
xy vanishes in the expansion of (6). This coefficient is found, after 
some trigonometric transformations, to be 


(x—@) sin (a—o.) + ycos (a ycosw,]sin (a+ f)| 


sin 2asin (8 —w,—a@,). 
Hence the horopter is symmetric when 
B— @, — @ = 0 (7) 
This is true, in particular, when the fixation is symmetric and 


f£ = 0. In general the locus of fixation points for symmetric horopters 
is the curve 


(2? — y? — a’) tan 6+ 2xy = 0, (8) 


which is a hyperbola whose center is at the origin and whose trans- 
verse axis has the inclination 1/4 — §/2. 

For symmetric convergence at the point [0, a tan (a —a,)], 
where w, = — w:=, the tangent to the horopter at this point is 


y = (x tan f cota —a) tan (a+), 


so that the slope here is tan f cot a tan (a—). All these tangents 
pass through the point (—a tan acot/, 0), which is the intercept of 
the linear horopter with the x-axis. 

Essentially three assumptions are involved in the derivation of 
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equation (6). One is the empirical fact that some horopter curve is 
linear. Another is that accommodative movements can be neglected. 
Presumably the validity of this assumption can be estimated by a 
consideration of the optical system of the eye. 

The third assumption is to the effect that the relation of corres- 
pondence is a fixed relation. This cannot be verified directly, but can 
only be inferred from theoretical considerations — unless, indeed, 
histological examinations can reveal the individual fiber pathways to 
the fusion center, wherever this may be. The theoretical neural mech- 
anism proposed by Verhoeff (1925) seems to imply such a fixity. And 
if equation (6) is found to agree with experiment, this fact would lend 
indirect empirical support to the notion. 

Actually equation (6) is based upon a further approximation, 
for it has been made on the assumption that the centers of the pencils 
of rays were also the centers of rotation of the two eyes. This, of 
course, is not the case, though the discrepancy would be negligible un- 
less the fixation point is near. The exact equation to replace (6) can 
be worked out in the same way, but is considerably more complicated. 

This work was aided by a grant from the Dr. Wallace C. and 
Clara A. Abbott Memorial Fund of the University of Chicago. 
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THE GENERAL FLUID CIRCUIT THEORY OF ACTIVE 
CHLORIDE ABSORPTION 


H. C. PETERS 
DEPARTMENT OF PHYSIOLOGY, UNIVERSITY OF TENNESSEE, MEMPHIS 


The original fluid circuit theory used to explain active intesti- 
nal absorption of chloride is modified to include diffusion and secre- 
tion of chloride and osmosis. The general differential equation de- 
veloped is integrated in a particular case. The definition, “effective 
concentration of chloride in the fluid passing into the intestinal 
lumen,” leads to simplified general expressions. 


In an attempt to explain the absorption of chloride against con- 
centration gradients which occurs in the ileum of the small intestine, 
R. C. Ingraham, H. C. Peters, and M. B. Visscher (1938, 1939) pro- 
posed the “fluid circuit” theory, which assumes that the absorption of 
water carries chloride out of the intestinal contents without changing 
its concentration from that in the lumen while a fluid free of chloride 
moves from the blood into the intestinal lumen. With the additional 
assumption that the rates of water transfer are constant the theory 
led to the equations: 

C/Cy= (V/Vo) *” (1) 
and 
V=V,—Dt, (2) 


where C and V are concentration and volume at time t, Co and Vp 
their initial values, R; the rate of movement of fluid into the intestinal 


- lumen, and D the rate of volume decrease. 


Concentration time curves for chloride absorption from a solution 
containing half isotonic sodium chloride and half isotonic sodium sul- 
fate agreed satisfactorily with the theory except at low chloride con- 
centrations. Since secretion and diffusion of chloride and osmosis have 


- been neglected in the derivation of equations (1) and (2), it seemed 


desirable to develop the theory on a more general basis. @his has been 
done as follows. We have 


I=CV, . (3) 
where I is the amount of chloride present in the intestinal lumen. 
dI/dt = C(dV/dt) + V(dC/dt) . (4) 
141 
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dI/dt =C,S + K(C, —C) —CR,, (5) 


where S is the rate of secretion of a fluid having a chloride concentra- 
tion C,, K a coefficient, C, the plasma chloride concentration, and Ry 
the rate of fluid movement out of the intestinal lumen. 

From equations (4) and (5) , 


C(dV/dt) + V(dC/dt) =C,S + K(C,—C) —CR,. (6) 


This is the equation for the general fluid circuit theory. It can be 
easily integrated, in the particular case in which dV/dt, C,, S, K, 
C,, and R, are constants, as follows: 

If dV/dt =—D, then, from equation (6), 


dC _CS + K(C,—C) —CR, + DC 


dt V ge 
Substituting dV/—D for dt and rearranging, 
ac a (8) 
(R, + K—D)C—CS—KC, DV~- 
Integrating between the limits C, and C, and V, and V, 
(R, + K — D)C — CS — KC, = ee (Ry+K-D)/D (9) 
(R, + K —D)C, —CS — KC, Vo 4 


where V is given by equation (2). 

Equations (9) and (2) will fit the data of Ingraham et al. satis- 
factorily at low chloride concentrations. At the present time however 
its use appears to be limited to relatively short experiments under 
restricted conditions. 

Equation (6) may be simplified by using the definition 


C.=[CS + K(C,;—C)]/R,, (10) 
where C, is the effective concentration of chloride in the fluid passing 
into the intestinal lumen at rate R;. R; includes osmotic as well as 
secretory water transfer. Substituting in equation (6), 

C(dV/dt) + V-(dC/dt) =CoR;— CRs. (11) 
At any instant 
dV /dt = R, — Ry, — Rez; (12) 


where R,,; isthe rate of water transfer by osmosis in the direction in- 
testine to blood only. If osmosis is in the opposite direction it is includ- 
ed in R; 7 

From equations (11) and (12), 
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V (dC/dt) =C. Ri —C(R; — Ry.) . (13) 
If dC/dt=0, then C,—C(R;—R,.)/R; and C.=C. 


Similarly if dC/dt < 0, then C. < C; and if dC/dt > 0, then 
Cs = C(R; ae, Ros) / Ri . 
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